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In  the  ionospheric  model  of  Venus,  the  observed  microwave  radiation  is  attributed 
to  free-free  emission  of  electrons  in  a  dense  Cytherean  ionosphere.  The  present  paper 
discusses  the  mechanisms  for  formation  of  such  a  dense  ionosphere,  both  in  the  original 
formulation  of  Jones  and  in  later  formulations  which  introduce  holes  in  the  ionosphere 
to  achieve  consistency  with  the  observed  radar  reflectivities.  Ionization  by  solar  ultra¬ 
violet  radiation  and  by  the  solar  proton  wind,  as  measured  near  Venus  by  Mariner  2,  are 
considered,  assuming  that  the  Cytherean  surface  magnetic  field  strength  is  <10~3  gauss. 
Both  sources  of  ionization  arc  primarily  effective  at  the  level  where  the  neutral  particle 
density  is  ^  10"  cm'3,  somewhat  above  the  level  at  which  diffusive  equilibrium  is  estal>- 
lished.  Three-body  eleetron-ion  recombination  processes  are  ineffective  at  these  densities. 
Radiative  recombination  is  plausible  only  if  all  the  X*  at  this  level  is  converted  to  N  by 
the  Herzberg-IIerzberg  photodissoeiation  mechanism.  But  the  times  for  X  to  diffuse  to 
dense  levels  where  recombination  to  X*  occurs  are  so  much  shorter  than  the  time  for 
Herzberg-IIerzberg  photodissoeiation  on  Venus,  that  X2  must  be  a  predominant  atmos¬ 
pheric  constituent  at  the  level  of  the  Cytherean  ionosphere.  If  dissociative  recombination 
therefore  prevails,  eharaeteristie  values  of  the  electron  density  are  n,  ~  10*  cm  3  over  a 
region  of  50-km  thickness.  If  radiative  recombination  prevailed,  w,  —  10s  ciir1  would 
result,  but  even  this  is  too  small  for  the  ionospheric  model  to  be  a  valid  explanation  of  the 
Venus  microwave  emission.  The  ionospheric  model  is  accordingly  rejected.  We  conclude 
that  the  observed  microwave  emission  arises  from  the  Cvthercaii  surface. 


Introduction 

An  extensive  set  of  observations  of  the 
brightness  temperature  of  Venus  as  a  func¬ 
tion  of  wavelength  and  phase  have  now  been 
performed  [see,  e.g.,  Mayer  (1903),  Barrett 
and  Staeliu  (1904),  Pollack  and  Sagan 
(1905)].  Near  inferior  conjunction,  tem¬ 
peratures  of  about  230°K  are  observed  in 
the  8-13  p  infrared,  rising  to  350°K  in  the 
4-8  mm  region,  and  then  to  about  580°K  at 
centimeter  and  decimeter  wavelengths.  A 
phase  effet  t  at  cent  imeter  wavelengths  of  the 
order  of  50  K°  now  appears  established,  ♦he 
dependence  on  phase  being  in  the  sense  of 
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an  increase  in  temperature  at  phase  angles 
removed  from  inferior  conjunction. 

Observations  of  the  radar  return  from 
Venus  at  wavelengths  of  08  cm,  43  cm,  and 
12.5  cm  by  Koteinikov  (1901),  Pettcngill 
cl  al.  (1902),  and  Victor  and  Stevens  (1901), 
respectively,  have  determined  that  the 
Cytherean  reflectivity  is  some  10-12%  the 
expected  return  from  a  perfectly  conducting 
sphere  of  the  same  dimensions.  Later  meas¬ 
urement  shave  confirmed  these  values.  There 
is  no  compelling  evidence  for  a  wavelength 
dependence  of  reflectivity,  and  the  radar 
results  have  Iteen  interpreted  as  indicating 
reflection  from  the  solid  surface  of  the  planet. 

Correlation  of  the  variations  from  the 


105 


©  1966  by  Academic  Press  Ine. 


V 


' 

-xta  _ 


106  R.  G.  WALKER  AND  C.  SAGAN 


mean  value  of  the  Astronomical  Unit  meas¬ 
ured  at  68  cm  by  Pettengill  el  al.  (1962)  with 
the  10.7-em  solar  flux  has  been  noted  by 
Whipple  (1962),  while  the  same  correlation 
with  the  20-cm  solar  flux  has  been  reported 
by  Priester,  Roemer,  and  Schmidt-Kaler 
(1962).  The  observed  correlation  is  in  the 
sense  that  an  increase  in  solar  flux  is  accom¬ 
panied  by  a  decrease  in  the  value  of  the 
Astronomical  Unit. 

Three  principal  models  have  been  pro¬ 
posed  to  explain  the  microwave  observa¬ 
tions.  The  greenhouse  model  (Sagan,  1960; 
Sagan  and  Pollack,  1966)  and  the  aeolosphere 
model  (Opik,  1961)  both  attribute  the  600°K 
radiation  to  emission  from  the  planetary 
surface.  In  the  ionospheric  model  (Jones, 
1961;  Sagan,  Siegel,  and  Jones,  1961),  the 
high  brightness  temperature  is  attributed  to 
free-free  transitions  of  electrons  in  a  deep 
Cytherean  ionosphere.  Depending  upon  the 
particular  model  chosen,  the  atmospheric 
structure  and  surface  conditions  on  Venus 
could  be  widely  different.  Critical  discussions 
of  the  relative  merits  of  these  models  have 
been  given  by  Kellogg  and  Sagan  (1961), 
Sagan  and  Kellogg  (1963),  and  Sagan  and 
Pollack  (1960).  In  the  preset  it  paper  we 
attempt  a  critical  appraisal  of  the  ionospheric 
model. 

Existing  Ionospheric  Models 

The  ionospheric  models  discussed  by  Jones 
(1961)  and  by  Sagan,  Siegel,  and  Jones 
(1961)  consider  the  free-free  emission  of 
electrons  in  an  isothermal  layer  at  7  e  = 
G0U°K.  If  I lw?  optical  depth  in  the  liver  at 
3  cm  is  about  2,  a  good  fit  to  the  microwave 
data  is  obtained.  The  dependence  of  the 
optical  depth  on  the  square  of  the  wave¬ 
length  renders  the  layer  optically  thick  at 
wavelengths  greater  than  5  cm,  and  trans¬ 
parent  for  wavelengths  less  than  1  cm.  To 
obtain  the  necessary  opacity  at  3  cm  we 
require  the  integral  of  the  square  of  the 
electron  density  through  the  layer  to  be  of 
the  order  of  4  X  1025  era-5  on  the  dark  side 
of  the  planet.  The  plkase  variations  reported 
at  8  mm  can  be  reproduced  if  the  electron 
density  is  two  to  five  times  greater  on  the 
sunlit  side.  A  major  difficulty  with  the 
ionospheric  model  is  the  maintenance  of  the 


high  electron  densities  required,  ne  ~  109  to 
1010  cm-1.  One  must  seek  mechanisms  for 
producing  a  very  high  ionization  rate  or  else 
invoke  extremely  low  values  of  the  electron- 
ion  recombination  coefficient.  Difficulties 
also  arise  in  explaining  the  radar  reflec¬ 
tivities  since,  if  reflection  arises  from  the 
surface,  the  high  opacity  of  the  ionospheric 
layer  would  cause  almost  complete  absorp¬ 
tion  of  the  signal.  To  overcome  this  difficulty 
Kellogg  and  Sagan  (1961)  have  suggested 
that  there  might  lie  a  “hole”  in  the  iono¬ 
sphere  near  the  antisolar  point  where  the 
electron  density  would  be  expected  to  decay 
during  the  long  Venus  night.  Priester, 
Roemer,  and  Schmidt-Kaler  (1962)  have 
suggested  that  the  radar  reflectivity  could 
be  explained  if  the  electrons  were  in  clouds 
of  A'e  ~  109  cm~3  covering  one-fourth  of  the 
Cytherean  surface  with  a  mean  density  of 
Ne  ~  10s  cm-3  in  the  regions  between  the 
clouds.  Thus  the  radar  return  would  lie  from 
the  surface  while  the  emission  would  arise 
from  the  electron  clouds.  The  electron  tem¬ 
perature  on  this  mix  lei  would  lie  800°K  for 
an  ionosphere  50  km  thick.  A  similar  sug¬ 
gestion  has  been  made  in  two  recent  papers 
by  Danilov  and  Yatscnku  (1963). 

In  their  first  paper  Danilov  and  Yatsenko 
conclude  that  the  ionospheric  model  can  lie 
reconciled  with  the  observations  if  radiative 
recombination  of  atomic  ions  is  the  dominant 
recombination  process  in  the  Venus  iono¬ 
sphere  However,  to  maintain  the  high 
electron  densities  observed  on  the  dark  side 
of  the  planet  an  unknown  and  enigmatic 
ionizing  agent  is  postulated.  No  evkkawe  is 
offered  to  support  its  presence.  It  is  assumed 
that  this  agent  is  similar  to  that  invoked 
to  account  for  the  observed  ionization  of  the 
Earth’s  F  layer  over  the  polar  regions  during 
the  long  polar  night.  The  radar  data  may  be 
explained  if  it  is  assumed  that  68-cm  reflec¬ 
tions  result  from  the  region  of  maximum 
electron  density  in  the  ionosphere,  43-cm 
radiation  is  totally  reflected  from  the  solid 
surface  and  attenuated  by  the  ionosphere, 
and  12.5-crn  radiation  is  partially  absorbed 
upon  reflection  from  the  planetary  surface. 

In  order  to  obtain  the  necessary  trans¬ 
parency  at  43  cm  while  still  retaining  the 
ionospheric  reflection  at  68  cm,  the  authors 
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require  the  electron  temperature  to  Ik.*  very 
high,  T,  ~  5000°K.  Such  a  higli-t  emperat  ure 
layer  would  exhibit  a  brightness  temperature 
rapidly  increasing  with  wavelength.  To  re¬ 
duce  this  abrupt  rise  and  thus  bring  the 
predicted  bright  ness  t  emperat  ures  int  o  rough 
agreement  with  those  observed,  an  overlying 
low  temperature  (Te  =  (XX)°Iv)  layer  is 
invoked.  The  ad  hoc  nature  of  these  assump¬ 
tions  is  quite  striking. 

In  their  second  paper  Danilov  and  Yat- 
senko  propose  a  mwlel  similar  to  that  of 
Priester,  Roemer,  and  Schmidt-Kaler  (1962). 
The  lower  ionosphere  (they  retain  the  two- 
laver  concept)  is  porous — there  are  clouds  of 
high  electron  density  and  interstitial  regions 
of  low  electron  density.  The  radar  signals 
pass  through  the  low-density  regions  un- 
attenuated  and  are  reflected  from  the  sur¬ 
face.  An  optimum  choice  of  parameters 
which  yield  a  reasonably  good  fit  to  the 
observed  spect  rum  is  a  temperature  of 
400°K  at  the  surface,  1500°K  in  the  hot 
layer,  and  700°K  in  the  cool  layer.  The 
surface  reflectivity  is  taken  to  be  0.4,  and 
the  clouds  of  high  electron  density  cover 
half  of  the  planetary  surface.  An  optical 
depth  of  2  at  X  —  10  cm  is  derived  for  both 
of  these  ionospheric  layers. 

Most  recently,  Kuzmin  (1964)  has  pro¬ 
posed  a  range  of  ionospheric  models  with 
holes.  He  explicitly  considers  the  effects  of 
collisions  with  neutral  particles  in  a  some¬ 
what  unusual  ionosphere,  where  the  collision 
frequency  greatly  exceeds  the  free-free 
emission  frequency  so  the  absorption  coeffi¬ 
cient  is  wavelength-independent.  The  ob¬ 
served  radar  reflectivity  will  not  necessarily 
be  the  true  reflectivity  of  the  surface,  be¬ 
cause  of  absorption  in  the  ionosphere  in  the 
interstitial  regions  between  the  holes.  If  the 
holes  cover  only  a  small  fraction  of  the  disk 
of  Venus,  the  true  surface  reflectivity  must 
be  very  large.  Thus,  for  a  smooth  surface 
with  20%  ionospheric  cover,  an  electron 
temperature  of  1300°K,  a  surface  tempera¬ 
ture  of  400°K,  and  a  dielectric  constant  of 
about  5  are  adequate.  On  the  other  hand,  if 
the  disk  is  70%  covered  by  an  ionosphere, 
the  electron  temperature  may  be  as  low  as 
650°K,  and  the  surface  temperature,  about 
500°K;  but  the  surface  dielectric  constant 
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must  now  be  about  15.  As  the  total  iono¬ 
spheric  cover  liecomes  greater,  the  required 
values  of  the  surface  temperature  rapidly 
exceed  those  required  if  the  microwave 
emission  arises  mostly  from  the  surface;  at 
the  same  time,  the  dielectric  constants  be¬ 
come  exceedingly  large. 

It  is  clearly  difficult  to  achieve  high  elec¬ 
tron  densities  in  a  region  of  high  neutral 
particle  densities.  As  we  shall  see  below, 
neither  solar  electromagnetic  radiation  nor 
solar  corpuscular  radiation  is  likely  to  pro¬ 
duce  high  electron  densities  in  regions  which 
have  neutral  particle  densities  in  the  range 
required  by  Kuzmin,  who  desires  1010  to  10" 
collisions/sec.  Such  a  collision  frequency  is 
characteristic  of  the  terrestrial  atmosphere 
near  the  surface.  The  corresponding  atmos¬ 
pheric  density  on  Venus  is  achieved  just 
below  the  clouds  (Sagan,  1962).  But  the 
Cytherean  atmospheric  rotational  tempera¬ 
ture  at  this  depth  is  300°K  or  less  (Sagan, 
1962),  not  the  >600°K  which  must  be  the 
electron  temperature  for  the  ionospheric 
model  to  function.  As  the  ionosphere  is 
moved  to  deeper  (and  hotter)  levels  the 
situation  is  improved,  but  the  whole  point 
of  the  ionospheric  model  is  to  avoid  a  hot 
lower  atmosphere.  With  no  other  source  of 
ionization  at  hand,  Kuzmin  invokes  either  a 
primary  cosmic-ray  flux  at  Venus  or  a  radio¬ 
active  particle  density  below  the  Venus 
clouds  several  orders  of  magnitude  larger 
than  the  corresponding  values  for  Earth. 
Those  models  which  place  the  ionosphere  in 
the  lower  Cytherean  atmosphere  are  too 
riddled  with  contradictions  and  improb¬ 
abilities  to  be  considered  further  here. 
Kuzmin’s  models  which  place  the  ionized 
layer  in  the  upper  Cytherean  atmosphere 
give  electron  densities  of  the  order  of  2  X  109 
cm-3,  in  agreement  with  previous  discussions. 

There  are  two  general  objections  to  the 
ionospheric  model  besides  the  ones  arising 
from  economy  of  hypothesis.  Drake  (1962) 
has  pointed  out  that  as  the  wavelength  of 
observation  decreases,  the  brightness  tem¬ 
perature  observed  in  an  optically  thick 
ionosphere  corresponds  to  greater  and 
greater  depths  within  the  ionosphere  Unless 
the  ionosphere  is  almost  exactly  isothermal — 
a  circumstance  of  small  probability — the 


108 


H.  (i.  WALKElt  AND  C.  SAC  AX 


centimeter  wavelength  brightness  tempera¬ 
ture  should  lx*  strongly  wavelength-de¬ 
pendent,  contrary  to  observation. 

Secondly,  an  isothermal  ionosphere  In¬ 
coming  optically  thin  at  a'oout  1  cm  should 
exhibit  limb-brigbte.iing  at  that  wavelength 
(Sagan,  Siegel,  and  Jones,  1901).  The 
Mariner  2  microwave  radiometer  channel  at 
19  mm  showed,  instead,  limb-darkening 
(Barath  el  al .,  1903).  Danilov  (1904)  has 
constructed  a  porous  ionospheric  model 
which  yields,  for  wavelength  •  near  1  cm. 
limlwlarkening  over  most  of  the  disk  and  a 
sharp,  thin  intensity  spike  at  the  limb.  While 
the  topographical  resolution  of  the  Mariner 
2  microwave  experiment  was  small,  it  never¬ 
theless  was  adequate  to  exclude  most  plausi¬ 
ble  ionospheric  spikes. 

But  there  seems  little  doubt  that  if  suffi¬ 
cient  ingenuity  were  employed,  very  complex 
models  could  be  invented  which  have  iono¬ 
spheric  hoios  strategically  placed  to  simulate 
limb-darkening,  and  at  the  same  time  have 
closely  isothermal  interstitial  material.  The 
key  property  of  all  ionospheric  models  is  the 
high  electron  density  required.  Since  this  is 
the  most  vulnerable  feature  of  such  models, 
it  appears  justified  to  examine  in  as  much 
detail  as  possible  what  conditions  exist  in 
atmospheric  regions  appropriate  to  the 
formation  of  a  Cythereau  ionosphere.  It  is 
the  purpose  of  this  paper  to  discuss  the 
mechanisms  by  which  an  ionosphere  may  Ik* 
formed  and  maintained  on  Venus,  with 
particular  emphasis  on  ionization  by  solar 
protons,  a  mechanism  first  suggested  by 
Jones  (1901). 

Basic  Relations 

The  specific  intensity  l{v,  8)  of  radiation 
emerging  at  angle  8  to  the  normal  in  a 
planetary  atmosphere  is  given  bv  the  formal 
solution  to  the  equation  of  radiative  transfer, 


H  =  cos  8.  and  the  optical  depth  t  is  given  by 


l(v,  8)  ■=  I  ( v,  7’,  )c~r " 


l,\v,  8, 


t)i-tIi‘iIt/h ,  (1) 


where  /,  is  the  specific  intensity  arising  from 
the  surface  at  emission  temperature  T„  Ia  is 
the  intensity  from  the  atmospheric  layers, 


-  /‘a 


( v ,  h)  dh, 


where  h  is  the  height  in  the  atmosphere,  and 
k{v,  h)  the  absorption  cross  section  per  unit 
volume.  At  radio  frequencies  the  Rayleigh- 
Jeans  approximation  to  the  blacklxxly  law 
may  lx*  used  to  express  the  intensity  as 

/(,,  T)  =  (2hir/c2)T,  (3) 

where  h  is  Planck’s  constant,  and  c  is  the* 
velocity  of  light.  Substitution  of  (3)  into  (1) 
and  integration  over  all  solid  angles  yields 
radio  brightness  temperatures, 


r„  =  2  f 

J  o 


T(ii)  m  dii, 


Tn  =  -T,Ei{rm)  +2  JQ  ^  T{r)e^drdn, 


where  r„,  refers  to  the  total  thickness  of  the 
at  mosphcric  layers.  Foi  t  tie  special  case  of  an 
isothermal  atmospheric  layer  aL  a  tempera¬ 
ture  7’„,  Eq.  (5)  reduces  to 

Tn  =  2T,i:,{rJ  +  Ta[l  -  2A’,(rm)J,  (0) 

where  E3(r)  is  the  third  exponential  integral, 
tabulated,  e.g.,  by  Kourganoff  (1952). 

The  absorption  coefficient  k  which  ap¬ 
pears  in  Eq.  (2)  is  given  by  Oster  (1901)  for 
free-free  transitions  of  electrons.  Adopting 
a  value  of  4.0  for  the  Gaunt  factor  yields, 
with  the  wavelength  X  in  centimeters, 

k  =  7.91  X  10-!3Xs(«e«,/77'2),  (7) 

where  n,  is  the  electron  density,  n,  the  ion 
density,  and  7’,  the  electron  temperature. 
The  optical  depth  from  Eq.  (2)  is  then 

r(X)  =  7.91 

fh 

X  I0-23X2  nr(h)n l(h)'I,r(h)~vi  dh,  (8) 
Jh  „ 

which,  for  the  special  case  of  an  isothermal 
layer  with  n,  =  «„  is 


r(X)  =  7.91  X  10~23X: 


27’r3/2  [  He 

j  ha 


2  dh.  (9) 
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Values  oe 

TABLE 

/  Iir  ilh  FOR  T 

1° 

=  1  AN1>  T, 

X  T.(fi) 

J it**  <ih  at  X  = 

T. 

0.5 

1.0 

3.4 

10 

21 

43 

fiS 

(MX) 

7.42 

1  SO 

1.01 

1.80 

4.21 

1  00 

4.01 

son 

114 

2.  SO 

2.48 

2.  SO 

0  50 

1 . 55 

0.18 

1000 

10. 0 

4. (HI 

3.40 

4(H) 

0.07 

2  16 

8.63 

1500 

29.3 

7.34 

0.30 

7  34 

10.7 

3.07 

15.8 

2000 

45.2 

11.3 

9.7S 

113 

25  7 

O  12 

24.4 

2500 

03.1 

15.8 

13.0 

15.8 

35 .  S 

S  54 

34.1 

3000 

S2.8 

20.7 

17.9 

20  7 

47  1 

11  2 

44  S 

Column 

X 

X 

X 

X 

X 

X 

X 

multiplier 

102‘ 

102< 

102* 

102* 

1023 

l()» 

1022 

a  /  n,2  ilh  in  cm-*;  X,  cm;  and  Te,  °K. 


Table  I  presents  values  of  the  integral  of  the 
square  of  the  electron  density  required  to 
achieve  unit  optical  depth  for  a  number  of 
selected  wavelengths  and  temperature  as 
calculated  from  Eq.  (9).  It  is  obvious  from 
the  table  that  an  increase  in  the  temperature 
of  the  layer  requires  a  substantial  increase 
in  the  integrated  electron  density  to  achieve 
the  same  optical  depth. 

Consider  now  an  atmosphere  in  hydro¬ 
static  equilibrium  satisfying  the  perfect  gas 
law.  The  pressure  P  at  any  level  is  deter¬ 
mined  by  the  relation 

dP/P  =  -  (mg/RT)  dh,  (10) 

where  g  is  the  local  acceleration  of  gravity, 
m  is  the  mean  molecular  weight  in  atomic 
mass  units,  T  is  the  temperature  in  °Iv,  and 
R  is  the  universal  gas  constant  [8.314  X 
107  erg  gm-1(K0)-1].  For  the  atmospherie 
models  considered  here  it  will  be  assumed 
that  g  and  m  remain  constant  for  all  heights 
h  above  the  height  hb  of  the  base  or  reference 
level. 

For  the  case  of  an  isothermal  atmosphere 
at  a  temperature  T,  Eq.  (10)  integrates  to 

P(h)  =  Pb  exp  [-(mg/RT)(h  -  hb)],  (11) 

where  Pb  is  the  pressure  at  the  reference  level 
hb. 

Integration  of  Eq.  (10)  for  a  temperature 
distribution  of  the  form 

T(h)  =  Tb  +  L(h  -  hb)  (12) 


and  L  ^  0,  yields  for  the  pressure 

n+rthr^J  •  <13> 

where  the  constant  L  is  the  temperature 
gradient. 

The  redueed  height  of  the  atmospheric 
layer  u’(/i)  is  obtained  from 

w(h)  =  22.4  X  10 z[P{h)/mg]  em-atm,  (14) 

where  P{h)  is  given  by  either  Eq.  (11)  or 
Eq.  (13)  depending  upon  the  model  chosen. 
If  we  are  interested  in  the  redueed  height 
for  the  tth  atmospheric  constituent,  then 
P(h)  is  replaced  by  the  partial  pressure 
P<(h). 

The  mean  molecular  weight  m  for  any 
given  atmospheric  layer  is  given  by 

m  =  ^  n,Wi/J  «,,  (15) 

t  * 

where  »,•  is  the  number  density  of  the  ith 
speeies  of  molecular  weight  m,-.  The  molecu¬ 
lar  weight  is  thus  determined  by  the  relative 
concentration  of  moleeular  and  atomic 
species  at  the  altitude  h,  and  thus  will  be 
dependent  upon  the  degree  of  dissociation  at 
that  altitude. 

Of  interest  is  the  frequency  with  whieh 
particles  of  different  types  eollide  at  any 
given  level  in  the  atmosphere.  The  collision 
frequency  v  is  determined  by  the  mean 
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velocities  r  of  the  particles  and  the  local 
mean  free  path  x,  through  the  relation 

v  =  ijx.  (10) 

For  a  Maxwell-Bolt zmann  distribution  of 
velocities 

v  =  (8  RT/rm)'13,  (17) 

and 

a-  =  (2  l«na»)-1;  (18) 

r  is  the  effective  collision  radius  and  is 
strongly  dependent  upon  the  types  of  collid¬ 
ing  particles.  For  collisions  between  atoms 
and  molecules  in  air  r  is  of  the  order  of  1.8  X 
10-8  cm  (Minzner  and  Ripley,  1957).  Sub¬ 
stitution  of  numerical  values  into  Eq.  (10) 
yields  for  collisions  in  air 

Vn  =  2.1  x  lO~u(T/m)'linH  sec'1,  (19) 

where  nH  is  the  total  number  density  of 
neutral  particles. 

Ginzburg  (1901)  has  discussed  various 
collisional  processes  appropriate  to  free  elec¬ 
trons  in  ionospheric  layers.  The  total  elec¬ 
tron  collisional  frequency  ve  is  given  by 

Ve  —  Ve-«  +  Ve-i  +  Ve-e,  (20) 

where  the  subscript  (e-n)  refers  to  collisions 
of  electrons  with  neutral  particles,  (e-i)  to 
collisions  with  positive  ions,  and  (e-e)  to 
electron-electron  collisions.  Substitution  of 
numerical  values  into  the  relations  derived 
bv  Ginzburg  yields  the  following  expressions: 

cc_„  =  6.3  X  10-»(7’/300)1/2n„  sec-1  (21) 

Vr-i  -  29.77’-:,,2«,[l  +  0.185  In  (T/ii"3)} 

sec-1  (22) 

v,.,  ~  ve-t.  (23) 

Equation  (22)  holds  for  regions  where  the 
electron  temperature  Te  is  equal  to  the  local 
kinetic  temperature  T.  For  regions  where 
Tc »  T,  Eq.  (22)  must  be  replaced  by  a  more 
elaborate  relation.  For  collisional  frequencies 
involving  positive  ions,  Ginzburg  derives 


wi-m  ^  1.5  X  10~10/f„( 77.300) 1,2  sec-'  (24) 

and 


where  mr  is  the  electron  mass  and  in,  the 
mass  of  the  positive  ion. 

For  some  purposes  it  will  lx?  convenient  to 
express  the  number  density  n (li )  in  terms  of 
w(h),  the  reduced  height  of  the  atmosphere. 
This  may  lx1  accomplished  by  combining 
Eq.  (11)  with  the  equation  of  state  giving 

n(h)  =  3.24  X  l0"mgT-'w(h).  (26) 

The  basic  equations  governing  the  time 
variation  of  the  electron  density  ne  in  an 
ionospheric  layer  are  discussed,  Tor  example, 
by  Mitra  (1952)  and  by  Xawroeki  and  Papa 
(1961).  These  relations  depend  upon  the 
rates  of  processes  which  tend  to  increase 
the  electron  density  (e.g.,  photoionization, 
photodetachment,  collisional  detachment, 
associative  ionization);  and  of  those  proc¬ 
esses  which  tend  to  decrease  the  electron 
density  (e.g.,  elect rou-ion  radiative  recom¬ 
bination,  dissociative  recombination,  elec¬ 
tron  attachment).  At  any  given  level  in  the 
atmosphere  it  is  usually  found  that  one  set 
of  processes  effectively  determines  the 
electron  density,  and  the  very  large  number 
of  other  contributing  reactions  can  be  neg¬ 
lected.  For  the  special  case  in  which  the 
electron  production  is  by  ionization  and  the 
decay  of  the  electron  density  is  controlled  by 
recombination  processes,  the  time  variation 
of  the  electron  density  is  given  by 

dnjdt  =  q  —  atne2,  (27) 

where  q  is  the  ionization  rate  (electrons 
cm-3  sec-1)  and  ar  is  the  effective  recombina¬ 
tion  coefficient  tem3  sec-1). 

The  decay  of  the  electron  density  once  the 
source  of  ionization  has  been  removed  is 
found  from  integration  of  Eq.  (27)  with 
q  =  0,  that  is 


nt(t) 


ne  o 

1  +  neoct'(t  -  to)’ 


(28) 


where  ne0  is  the  value  of  the  electron  density 
at  t  -  t0.  The  time  interval  At  required  for 
the  electron  density  to  fall  to  1/2  of  its  initial 
value,  obtained  from  Eq.  (28)  by  jetting 
ne(l)  =  (l/2)ne0,  is 

At  =  (n«a«)-'.  (29) 


The  quantity  At  is  thus  the  mean  lifetime 
of  an  electron. 
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At  equilibrium  in  the  layer,  dn,/dt  =  0, 
and  (27)  reduces  to 

q  =  a,n,-  (30) 

or 

*«  =  (7/«,)1/2.  (31) 

Table  II  presents  some  values  of  q  required 
by  Eq.  (30)  to  satisfy  various  combinations 
of  ae  and  nc. 

TABLE  It 


Valves  of  rm:  Ionization  Kate  q  Required  tsv 
Eg.  (30) 


a* 

(cm1  sec"1 

V  (cm’ 

)  at  «c  (cm 

-1)  - 

)  I01 

I0« 

10" 

10' 

10* 

!0» 

to-* 

10-' 

10* 

10‘ 

10s 

1010 

10” 

10-10 

1 

10- 

10* 

10‘ 

108 

10‘» 

10-‘J 

10- 

1 

105 

10< 

10‘ 

10* 

10-“ 

10-« 

10-! 

1 

101 

104 

10‘ 

Photochemistry  of  the  Upper 
Atmospheric  Layeils 

It  is  of  interest  to  consider  briefly  what 
may  lie  deduced  concerning  conditions  in 
the  upper  atmosphere  of  Venus.  Observa¬ 
tions  of  the  occultation  of  Regulus  by  de 
Vaucouleurs  and  Menzel  (1960)  indicate  a 
pressure  at  the  occultation  level  of  2.6  dynes 
cm-2.  These  observations  together  with 
several  observations  pertaining  to  the  tem¬ 
perature  and  scale  height  at  this  level  have 
lieen  discussed  by  Sagan  (1962),  who  tenta¬ 
tively  concluded  that  the  following  condi¬ 
tions  exist  at  this  level:  To  =  203°K,  mo 
=  28.8,  and  the  scale  height  H0  =  6.8  km. 
The  value  of  the  mean  molecular  weight  m0 
is  based  upon  a  mixing  ratio  of  a  C02  and 
(l  —  a)  X2  below  the  occultation  level, 
where  a  ^  0.05  and  the  C02  is  50%  photo- 
dissociated  at  the  occultation  level.  If  the 
C02  abundance  is  smaller  (Chamberlain, 
1965)  m0  will  decline  very  slightly.  Sub¬ 
stituting  the  above  values  into  Eq.  (14)  and 
adopting  g  =  870  cm  sec-2  we  find  the  re¬ 
duced  height  of  the  atmosphere  above  the 
occultation  level  to  be  u>n(N2  +  C02)  = 
2.32  em-atm,  U’U(N2)  =  2.21  cm-atm,  and 
w0(CO2)  =  0.11  cm-atm.  From  Eq.  (26)  the 
total  number  density  is  9  X  1013  cm-3,  and 


from  Eq.  (19)  the  collision  frequency  vn  = 
5  X  10*  sec-1. 

In  the  region  above  the  occultation  level 
photoionization  and  photodissociation  will 
be  important  proeesse®  determining  the 
upper  atmospheric  composition.  For  mono¬ 
chromatic  radiation  of  frequency  v  in  a 
spectral  region  where  the  absorption  by  only 
one  atmospheric  constituent  X  is  important, 
the  optical  depth  in  the  zenith  is  given  by 

r,(X)  =  k,(X)w(X),  (32) 

where  k,(X)  is  the  absorption  coefficient  in 
cm-1  arid  w{X)  is  the  reduced  thickness  of 
constituent  X.  In  terms  of  the  absorption 
cross  section  <r,(X),  k,(X)  =  N0<t,(X)  where 
No  is  Losehmidt’s  number,  2.6875  X  1019 
cm-3.  Substitution  into  Eq.  (32)  yields  the 
reduced  thickness  required  to  achieve  unit 
optical  depth, 

in, (X)  =  3.72  X  lO-a’MX)]-1  cm-atm.  (33) 

Absorption  cross  sections  for  important 
atmospheric  gases  have  been  published  by 
Nawrocki  and  Papa  (1961),  which  is  the 
source  for  all  uncited  cross  sections  and 
reaction  rates  in  the  following  discussion. 

Sagan  (1961,  1962)  has  pointed  out  that 
C02  is  dissociated  on  Venus  by  solar  ultra¬ 
violet  radiation  according  to 

C02  +  hv—*  CO  +  O  (\  <  1692  A),  (a) 

where  unit  optical  depth  occurs  near  the 
occultation  level  (i0i(CO2)  =  0.12  cm-atm 
for  r  =  3.0  X  10-19  cm2].  This  reaction  pro¬ 
ceeds  if  there  is  less  than  0.03  cm-atm  of 
02  above  the  occultation  level,  since  the 
Schumann-Runge  dissociation  continuum 
reaches  unit  optical  depth  for  this  thickness 
at  X  1692. 


Some  possible  reactions  involving  the 


products  of  (a)  are 

CO  +  O  +  Y  -»  C02  +  Y 

(b) 

O  +  O  +  Y— »Oi  +  Y 

(c) 

O  +  02  +  \  — >  O3  -(-  A 

(d) 

O  -)-  O  — *  02  -(-  hv 

(e) 

02  +  hv-*0  +  O(X1250  -  1760  A)  (f) 
o(3P)  +  co  — » eo2.  (g) 
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Reactions  (b),  (c),  and  (d)  are  three-body 
processes,  while  (e)  is  a  two-body  radia¬ 
tive  process.  Following  Mitra  (1932),  we 
find  that  if  N  is  the  total  particle  num¬ 
ber  density,  K  the  fraction  of  N  which 
is  the  principle  reactant  (e.g.,  O  in  the 
above  reactions),  and  (1  —  K) X  the  con¬ 
centration  of  the  third  body,  then  the 
rate  of  recombination  by  three-body  proc¬ 
esses  <7(3-body)  =  A2(l  —  K)X30,  where 
0  is  the  three-body  raie  coefficient.  The 
recombination  rate  by  two-body  processes 
is  </(2-body)  =  A'2A'2a,  where  a  is  the  rate 
coefficient  for  two-body  recombination.  Thus 
one  may  write 

q (3-body) /g(2-l tody)  =  (1  -  A) A' (d/a). 

(34) 

If  we  choose  K  «  1,  which  favors  three- 
bod}7  processes,  then  from  Eq.  (34)  it  is 
seen  that  for  A*  »  a/0  three-lrodv  processes 
will  dominate  and  conversely  for  .V  «  a/0 
two-body  recombination  will  dominate.  Rate 
coefficients  for  (c)  and  (d)  given  by  Xawroeki 
and  Papa  (19(il)  are 

0(c)  =  3  X  10-”  cm6  sec"1 

d(d)  =  2  X  10-34  cm6  sec-1. 

The  rate  coefficient  for  (e)  has  lx*en  esti¬ 
mated  as  a(e)  =  10_lfi  cm3  sec-1  if  both  0 
atoms  are  in  the  */*  state;  if  one  is  %P  and 
the  other  ll)  the  coefficient  may  he  as  high 
as  a(e)  =  10-14  cm3  see-1.  Coefficients  for 
(b)  have  not  lteen  found;  however,  it  is 
unlikely  they  will  differ  greatly  from  10-33 
to  10-34  cm6'  sec-1.  For  0  =  10-33  cm6  sec-1 
and  a  =  10-16  cm3  sec-1  we  have  Ar  »  1017 
cm-3  as  the  condition  for  three-bodv  proc¬ 
esses  to  dominate.  Since  A*  ~  1014  cm-3  at 
the  occultational  level  we  may  conclude  that 
three-body  collisions  are  unimportant  above 
the  occultation  level. 

The  rate  q  of  phot  odissociat  ion  is  deter¬ 
mined  by 

q  =  J  fjr,n  dv  (33) 

where  j,  is  the  solar  photon  flux,  <r,  is  the 
eross  section;  and  n  is  the  number  density. 
For  reaction  (a)  the  soiar  photon  flux  at 


A  <  1092  A  at  Venus  is  4  X  1012  photons 
cm-2  sec-1,  f^3X  I0~ ,s  cm2,  and  n(CO*) 
=  4.4  X  1012  cm-3.  Thus  the  rate  at  which  () 
is  formed  is  </(())  =  3.3  X  10®  atoms  cnr3 
sec-1.  If  (e)  proceeds  rapidly  the  equilibrium 
concentration  of  ()2  will  be  determined  by 
the  reactions  (a)  and  (f).  The  rate  of  de  ple¬ 
tion  of  Ojo  is  determined  by  absorption  of 
A  <  1750  A  in  the  Sehumann-Runge  con¬ 
tinuum.  Foi  this  reaction  the  flux  is  4  X  1012 
photons  cm2  sec-1,  a  -  1.8  X  10_IT  cm2,  and 
the  nite  q{()-)  -  7.2  X  I0~s  «(Os).  In  equi¬ 
librium  q{ O)  =  q{ Oj)  or  «(Os)  =  7  X  101" 
cm-3;  the  corresponding  reduced  thiekness 
for  Os  is  1.7  X  10~3  cm-atm.  This  represents 
an  upper  limit  to  the  0;  concentration  above 
the  occult ut ion  level.  If  (e)  proceeds  slowly 
the  amount  of  02  will  l>e  less  and  a  small 
equilibrium  abundance  of  ()  will  l>e  main¬ 
tained.  If  the  assumptions  inherent  in  the 
above  analysis  are  correct,  one  expects  to 
find  no  Schnmaim-Runge  shielding  of  COs 
and  no  appreciable  amount  of  ozone  above 
the  occultation  level.  This  analysis  holds 
only  for  the  sunlit  side  of  the  planet.  During 
the  long  Venus  ’light  photodissociation  ceases 
and  decay  of  CO  will  proceed  by  (b)  and  (g). 
while  the  abundance  of  O.  will  probably  bo 
determined  by  a  combination  of  (c),  (d),  and 
(e). 

It  is  to  be  noted  that  solar  radiation  at 
° 

A  <  1020  A  can  photoionize  Os,  the  maxi¬ 
mum  cross  section  being  a  =  4  X  10-isem2at 
A900.  The  reduced  height  for  optical  depth 
unity  is  about  10-2  cm-atm.  Thus  ionization 
can  occur  at  the  occulation  level.  Ionization 

o 

of  ()  can  occur  for  X  <  91 1  A  with  r  =  1  at 
m-’i(O)  =  10~2  cm-atm. 

Solar  radiation  in  the  spectral  region 
790  —  200  A  will  photoionize  X _>.  The  mean 
value  of  the  cross  section  for  this  region  is 
2  X  10-7  cm2.  At  unit  optical  depth  Un¬ 
reduced  thickness  is  ?C|(Ns)  =  2  X  10-3 
cm-atm.  Due  to  the  assumed  great  abun¬ 
dance  of  N*,  radiation  in  these  wavelengths 
will  not  penetrate  very  deeply  in  the  Cy- 
therean  atmosphere.  Other  absorptions  occur 
in  theoLyman-Birgc-IIopfiold  bands  (1450  — 
1000  A).  Cross  sections  for  the  strongest 
bands  are  on  t  he  order  of  10-20  cm2,  yielding 
tCi(L-B-H,  N2)  ~4  cm-atm.  Thus  radiation 
of  A  >  800  A  is  essentially  unattenuated  by 
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N.-  in  the  up|KT  atmosphere  alxne  tlie 
nmiltntion  level. 

Dissociation  of  N s  can  occur  hv  photo¬ 
ionization  and  dissociative  recombination, 

X*  +  hi*  -*  X  .+  +  e  (X  <  7%  A)  (h) 
X«+  +  c  -  X  +  X  (i) 

or  through  the  IlerzlxTg-IIerzlxM'g  predis¬ 
sociation  mechanism  (see  Bates,  1954), 

N*(X'2/)  +  hv  ->  Xs(a'lU  (j) 
(at  1200  —  1250  A)  followed  by 

X2(a'nff)  -*  X  +  N.  (k) 

In  the  abscnce'of  other  reactants  the  atomic 
nitrogen  recombines  by  a  t  hree-body  process, 

N  +  N  +  Y-»N,  +  Y.  (1) 

On  Venus  the  rate  for  reaction  (h)  is  (/(h)  ~ 
2.8  X  10~7  n(X«)  while  that  for  (j)  would 
lie  q(  j)  ^4X  I  O'2  **(Xa).  At  high  altitudes 
where  X790  A  can  penetrate  and  the  number 
density  is  low,  (h)  will  dominate,  while  at 
lower  altitudes  (j)  will  become  important. 
Due  to  the  slow  rate  of  (1),  /3(1)  ~  I0-27 
cm6  sec-1,  it  is  likely  that  nitrogen-oxygen 
lead  ions  may  produce  small  concentrations 
of  XO  and  X()». 

From  recent  balloon  observations,  Bot- 
tema,  Plummer,  and  Strong  (BHD)  estimate 
that  the  amount  of  water  vapor  above  the 
cloud  level  is  —  1 0  2  gm  cm-2.  Water  vapor 
has  a  weak  dissociation  continuum  at  A  < 

o 

1850  A.  The  maximum  cross  section  is  at 
A1050  where  (c, (lb/))  —  7  X  10-3  gm  cm-2. 
Thus  it  is  likely  that  lf»0  is  partially  dis¬ 
sociated  above  the  cloud  level,  in  the  absence 
of  protection  bv  ()*  or  other  absorbers. 

Ionization  Rate  Due  to  the  Solar 
Proton  Wind 

Jones  (1901)  has  suggested  that  solar 
protons  might  provide  a  significant  source 
oi  ionization  on  Venus.  This  suggestion  will 
be  considered  in  some  detail  here  for  an 
atmosphere  having  the  required  charac¬ 
teristics  of  the  ionospheric  model.  In  the 
following  discussion  it  is  assumed  that  Venus 
has  a  negligible  magnetic  field  strength,  an 
assumption  consistent  with  but  not  proved 
by  Mariner  2  magnetometer  results. 


The  depth  of  penetration  of  protons  into 
atmospheric  gases  lias  been  studied  by  Cook, 
Jones,  and  Jorgensen  (1953).  For  proton 
energies  less  than  70.5  k  *V  the  experimental 
data  fit  the  empirical  k  at  ion 

R  =  3.26E0  73,  (36) 

where  E  is  the  proton  energy  in  keV  and  R 
is  the  range  (in  cm)  of  the  protons  in  air 
at  1-mm  pressure  and  288°K.  The  reduced 
thickness  of  penetration,  ic,  in  cm-atm  (at 
STP)  is  then 

w  =  4.06  X  10~3E” 73  cm-atm.  (37) 

Values  of  the  proton  flux  and  energy 
spectrum  in  interplanetary  space  as  meas¬ 
ured  by  Mariner  2  have  been  reported  by 
Xeugebauer  and  Snyder  (19(52).  The  per- 
t  incut  data  from  their  report  are  summarized 
in  Table  III  together  with  the  range  w 


TABLE  III 

Ranoe  and  Eneroy  of  Interplanetary  Protons 


EQ 

(volts) 

r 

(km  st*c ) 

K 

(k.-V) 

IC 

(cm-atm) 

h' 

<%) 

510 

:tl4 

0.517 

0.0325 

— 

751 

:tst) 

0.700 

O.OU33 

18.3 

1124 

405 

1.14 

0.0045 

22.5 

1004 

50T 

1  05 

0  005*1 

30.5 

247(5 

O'.M) 

2.50 

0.(HI7!» 

10.0 

:ms 

S40 

3.70 

0  0100 

0.3 

540S 

1010 

5.42 

0.0130 

— 

S224 

1250 

S .  22 

0  0100 

— 

calculated  from  Eq.  (37)  and  the  proton 
energy  E  calculated  on  the  assumption  that 
all  particles  measured  were  protons.  The 
column  headed  F  gives  the  fraction  of  the 
time  the  maximum  of  the  spectral  flux  oc¬ 
curred  at  this  value  of  the  encrgy-to-charge 
ratio,  E/Q.  A  typical  value  for  the  particle 
density  during  periods  of  no  geomagnetic 
storms  was  given  as  2.5  protons  cm-3  at 
E/Q  =  1124  volts.  For  the  purposes  of  the 
following  discussion  the  values  of  the  proton 
flux,  WpP,  in  Table  IV  will  be  adopted  as 
representative  of  average  values  at  the  top 
of  the  Cvtherean  atmosphere.  This  spectrum 
assumes  np  =  4.0  protons  cm-3  for  751  volts 
<  E/Q  <  2476  volts,  np  =  0.4  cm-3  for 
E/Q  =  3688  volts,  and  np  =  0  outside  these 
limits.  These  values  of  np  were  obtained  from 
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the  Mariner  2  data,  and  adjusted  to  the 
Venus-Sun  distance  by  direct  scaling  of  the 
solar  wind  model  of  Brandt  (19(32).  It  is  be¬ 
lieved  that  the  adopted  proton  spectrum 
represents  a  reasonable  upper  limit  to  the 
average  proton  flux  at  Venus  during  quiet 
periods. 

TABLE  IV 
Solar  Proton  Flux 

E  (keV)  0.76  1.14  1.65  2.50  3.70 

npv  ( X  10’ cm-2  1.6  2.1  2.3  2.80  0.34 

sec-1) 


monoenergetie  beam  of  protons.  The  total 
ionization  rate  q(h)  may  then  be  found  by 
summing  over  the  proton  energy  spectrum. 
Adopting  a  value  of  30  eV/ionization  for 
Ei  and  inserting  numerical  values  for  the 
constants,  we  find  for  the  isothermal  layer 

g(A)  =  3.02^M^ 

) 

XeXp[“5?(*“A*)],  (43) 
and,  for  a  linear  temperature  gradient, 


The  ionization  rate  q(h,  E)  is  proportional 
to  the  rate  of  energy  loss  dE/dh  along  the 
path  of  the  proton,  and  to  the  proton  flux 
npv,  it  is  inversely  proportional  to  the  energy 
Ei  expended  per  ionization.  Accordingly, 

q(h,  E)  =  X  electrons  cm-3  sec-1, 

’  Ei  dh  10s 

(38) 


for  Ei  in  electron  volts  and  h  in  kilometers. 
The  rate  of  energy  loss  is 


dE  dE  dw 
dh  dw  X  dh 


(39) 


Differentiation  of  (37)  gives 

dE/dw  =  3.37  X  lO2/?0-27  kcV/cm-atm.  (40) 


#)  =  3.02  ^^(WptW-27 

ImglRL'l+l 

,  (44) 

where  the  summations  may  be  carried  out 
for  all  values  h  >  hj,„in  the  minimum  alti¬ 
tude  to  which  protons  of  energy  E,  can 
penetrate.  The  minimum  altitude  of  penetra¬ 
tion,  obtained  by  equating  Eq.  (14)  and  Eq. 
(37),  substituting  the  proper  expression  for 
P(h),  and  solving  for  h  -  hj,min  is 


T„ 
+  L(h 


Th  [  / (3.34  X  10»M*W- 
L  [\  mEj°-n  ) 


(45) 


For  an  isothermal  layer  dw/dh  is  found  by 
substitution  of  (11)  into  (14)  and  differen¬ 
tiating, 


dw 

Th 


22.4  X  10* 
RT 


Wt  <*  -«]• 

(41) 


For  a  layer  with  a  linear  temperature 
gradient  L  7*  0,  dw/dh  is  found  from  (13) 
and  (14), 


dw 

Th 


22.4  X  10* 

RTh 


X/J _ *± _ T 

A  *  .  Tk  +  L{h  -  A*)J 


mglRL)+\ 


(42) 


Multiplication  of  Eq.  (40)  by  Eq.  (41)  or 
Eq.  (42),  depending  upon  the  model  atmos¬ 
phere  being  considered,  and  substitution 
into  Eq.  (38)  yields  the  ionization  rate  for  a 


for  the  linear  temperature  gradient,  and 


hj,min  hb 


RT  / 6.34  X  103/V\ 
mg  \  mE/)n  )’ 


(46) 


for  the  isothermal  layer. 

Model  atmospheres  based  upon  the  above 
equations  have  been  constructed  to  examine 
the  rate  of  electron  production  on  Venus  by 
the  impact  of  solar  protons.  The  models 
presented  here  have  been  developed  subject 
to  the  constraint  that  most  of  the  proton 
energy  is  depleted  in  ionizing  the  particular 
atmospheric  region  of  interest,  a  procedure 
which  maximizes  the  derived  local  electron 
densities.  In  the  case  of  the  isothermal 
models  this  constraint  determines  the  pres¬ 
sure  Pb  at  the  base  of  the  isothermal  layer. 
From  Table  III  it  is  seen  that  3.7  keV 
protons  penetrate  to  0.0106  cm-atm.  Sub- 
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siitution  of  this  value  into  Eq.  (14)  yields  the 
base  pressure 

Pb  =  P{hm,n)  =  4.75  X  10~7  mg  dynes  era-2. 

(47) 

In  the  ease  of  a  linear  temperature  gradient 
the  above  const  raint  was  used  to  determine 
the  temperature'  gradient  L.  If  it  is  desired  to 
produce  electrons  at  a  mean  temperature 
ft  then  the  majority  of  the  protons  should 
be  thermalized  near  the  level  in  the  atmos¬ 
phere  where  T  =  Tt.  (In  this  and  the  above 
discussion  it  is  assumed  that  the  electrons 
are  quickly  thermalized.  This  assumption 
will  be  discussed  in  some  detail  later.)  From 
Table  III  a  weighted  mean  value  of  proton 
penetration  is  0.000  em-atm.  Substitution 
of  this  and  (15)  into  (14)  gives 

P  r  rp  ]<»»/KL) 

2.08  X  10~T  =  — '  ,  ~7 - j-r 

mg  [  /;,  +  L(h  —  hb)  J 

(48) 

By  taking  Tb  and  Pb  to  lie  the  temperature 
and  pressure  at  the  occultation  level  and  by 
substitution  of  Tt  foi  Tb  +  L(h  —  hk),  Eq. 
(48)  may  be  solved  for  L. 

Figures  1  through  4  show  the  ionization 


curves  obtained  for  several  model  atmos¬ 
pheres.  The  constant  parameters  for  the 
models  with  a  temperature  gradient  are 
To  =  203°K,  P<>  =  2.6  dynes  cm-1,  while  the 
variables  are  T  and  m.  For  the  isothermal 
models  the  ratio  T/m  is  the  variable  param¬ 
eter.  It  is  obvious  from  these  figures  that 
increasing  the  temperature  increases  the 
thickness  of  the  layer  of  electron  formation 
and  decreases  the  maximum  ionization  rate. 
The  jagged  or  multilayered  appearance  of 
the  curves  results  from  treating  the  Mariner 
2  flux  data  as  five  discrete  monoenergetic 
proton  fluxes.  In  reality  the  proton  spec¬ 
trometer  accepted  a  fairly  broad  energy  band 
AE  for  each  value  of  E/Q  measured.  It  is 
expected  that  proper  weighting  of  the  ob¬ 
served  flux  by  the  effective  slit  function  and 
smoothing  of  the  spectrum  to  account  for 
those  energies  falling  between  the  array  of 
fixed  slits  (and  not  subsequently  observed) 
would  lead  to  a  smooth  ionization  rate  falling 
within  the  envelope  defined  by  the  extreme 
points  in  the  figures. 

Defining  Q  as  the  integral  of  the  ionization 
rate  over  altitude 

Q  =  j  q(h)  dh .  (49) 


Fid.  1.  Ionization  rate  profile  due  to  solar  protons  measured  by  Mariner  2,  for  an  isothermal  layer  with 
temperature,  T,  and  mean  molecular  mass  m,  as  shown. 
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Fig.  2.  Ionization  rate  profile  due  to  solar  protons 
measured  by  Mariner  2,  for  a  layer  with  linear  tem- 
I>eraturc  gradient  L,  mean  electron  teiu|>eruturo 
t„  and  mean  molecular  mass  m,  as  shown. 


Fig.  3.  Ionization  rate  profile  due  to  solar  protons 
measured  by  Mariner  2,  for  a  layer  with  linear  tem¬ 
perature  gradient  L,  mean  electron  temperature  T,, 
and  mean  molecular  mass  m,  as  shown. 


Fig.  4.  ionization  rate  profile  due  to  solar  protons 
measured  by  Mariner  2,  for  a  layer  with  linear  tem¬ 
perature  gradient  L,  mean  electron  temperature 
7’,.  and  mean  molecular  mass  m,  as  shown. 


we  find  from  (30)  that,  :it  equilibrium, 

j  q(h)flh  =  j  a,.{h)Hr(h)2  <lh,  (50) 
which,  for  a,(h)  =  a,  =  constant,  gives 


Q  =  dr  j  nc(h)- rill.  (51) 

(J  may  lie  determined  from  integration  of 
Kq.  (38)  and  summation  over  the  proton  flux 


♦-sejo©.©* 

J 

Q  =  4.0  x  105  ^  (n„r), 


(52 1 
(531 


a  quantity  which  is  independent  of  the 
at mospherie  model.  .Substitution  of  the 
values  of  the  adopted  proton  flux  into  Eq. 
(53)  yields  Q  =  7.12  X  10,u  cm--’  sec  \ 
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A  quantity  useful  for  calculating  the  total 
opacity  of  the  electron  layer  is 

E  =  j  q{h)Tt{h)~2  -  dh,  (54) 

which  for  an  ionized  layer  in  equilibrium,  by 
Kq.  (30),  is 

S=  j  at(h)Tt(h)-3lint(hydh,  (55) 
and  for  at{h)  =  at  =  constant, 

E  =  a,  /  T'(h)-3lint(h)-  dh.  (56) 

For  the  isothermal  model 

E  =  QTr'n  =  7.12  X  101(Tr3/i  cm-  2  sec"1, 

(57) 

and  for  the  linear  temperature  gradient 
integration  of  Eq.  (38)  after  dividing  by 
Tt(h)m  yields 

1.37  X  103 

4  lURWmg)  +  1]  EiTf* 

1.81  X  10_7mgf 

1\ 

(58) 

As  it  pertains  to  the  total  opacity,  a  mean 
or  effective  electron  temperature  Tt  may  be 
defined  by 

frw  j  q(h)  dh  =  j  T,(h)~3liq(h)  dh,  (59) 

which  is  useful  where  ar(h)  ar. 

Values  of  the  chosen  parameters  m,  T„ 
and  the  derived  quantities  L,  T„  and  E  are 
presented  in  Table  V,  for  the  linear  models 
along  with  the  maximum  ionization  rate 
qm.  The  empirical  relation 

qm  =  5.5  X  \Qhn/Tr  cm-3  sec-1  (60) 

derived  from  the  table  holds  to  within  a 
maximum  error  of  7%. 

Similarly,  values  characteristic  of  the 
isothermal  models  are  presented  in  Table 
\T.  The  height  of  the  base  of  the  isothermal 
region  above  the  oceultation  level,  fo  —  ho, 


TABLE  V 

Characteristics  of  Model  Ionospheres  with  a 
Linear  Temperature  Gradient 


m 

i 

(K'Ani) 

T. 

UK) 

f. 

UK) 

qm 

(cm"1  sec-1) 

(cm-5  sec-1) 

28 

6.54 

600 

688 

2.34  X  104 

3.94  X  10* 

28 

4.76 

1000 

1210 

1.29 

1.69 

28 

3.90 

1500 

1890 

0.88 

0.87 

28 

3  62 

2000 

2550 

0.60 

0.55 

10 

7.04 

600 

686 

1.19 

3.95 

16 

5.10 

1000 

1182 

0.76 

1.74 

16 

4.26 

1500 

1860 

0.47 

0.89 

16 

3.86 

2000 

2520 

0.36 

0.55 

35 

6.32 

600 

682 

2.72 

3.98 

was  calculated  assuming  a  linear  tempera¬ 
ture  gradient  from  the  oceultation  level  to 
the  isothermal  region.  The  values  of  Table 
VI  fit  the  empirical  relation 

qm  =  4.7  X  IQhn/f,  cm-3  sec-1.  (61) 


TABLE  VI 

Characteristics  of  the  Isothermal  Models 


T. 

H 

Qm 

Z  hk  —  h 

UK) 

m 

(km) 

(cm-*  «ec_I) 

(cm-1  sec-*) 

(km) 

600 

2S 

20 . 5 

2.17  X  104 

4.84  X  10* 

68 

1050 

28 

37.5 

1.25 

2.09 

96 

600 

16 

37.5 

1.25 

4.84 

131 

Photoioxization  ox  Venus 

For  purposes  of  comparison  with  the 
ionization  due  to  solar  protons;  it  is  of  in¬ 
terest  to  obtain  an  order-of-magnitude 
estimate  of  the  ionization  rate  due  to  solar 
ultraviolet  radiation.  The  rate  of  photo¬ 
ionization  qp  for  radiation  in  the  zenith  is 
given  by 

qp  =  J  F /r,n  dv,  (62) 

where  F,  is  the  ultraviolet  flux  at  frequency 
v,  a,  is  the  cross  section  for  photoionization, 
and  n  is  the  number  density  of  particles  of 
cross  section  a,.  The  maximum  rate  of  ioniza¬ 
tion  will  occur  at  optical  depth  unity  in  an 
isothermal  layer  of  depth  given  by  Eq.  (33). 
For  monochromatic  radiation  at  the  top  of 
the  atmosphere  of  flux  V 

F,  =  F, 


^3 RLUmg 


r  '■  ■ 

vi 


(63) 
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The  maximum  rate  of  ionization,  qm,  is 

qm  =  <r,(X)n(X)F.ac-1,  (04) 

which,  by  virtue  of  Eqs.  (2(5)  and  (33),  with 
g  =  870  em  sec-2.  Incomes 

qm  =  3.9  X  10-<[wi(X)/7’]/,',o  cm-3  sec-1. 

(95) 

Solution  of  Eq.  (05)  with  T  =  1000°K  for 
X,  X2,  O,  and  ()»  yields  qm  =  7<>00,  15  000, 
0200,  and  25  000  electrons  cm-3  sec-1  at 
u>,  =  0.0034,  0.0015,  0.011,  and  0.0071 
cm-atm,  respectively.  Depending  somewhat 
upon  the  degree  of  dissociation  of  X2  the 
ionization  maxima  for  X  and  X»  will  fall 
within  the  region  of  significant  proton  ioniza¬ 
tion  (see  Table  HI).  Due  to  the  probable 
small  abundance*  of  O  and  ()2  the  maxima 
for  these  ionization  rates  will  most  likely 
fall  near  or  below  the  occultation  level.  [The 
values  for  the  ionization  cross  sections  and 
the  solar  ultraviolet  flux  were  taken  from 
Xawrocki  and  Papa  (1901).  The  solar  flux 
was  increased  a  factor  of  2  for  Venus.) 

A  generous  upper  limit  for  the  number  of 
electrons  formed  per  cm2  column  may  lx* 
obtained  by  assuming  tliat  all  photons  effec¬ 
tive  in  producing  ions  each  produce  one 
electron;  that  is, 

/  qP(h)  dh  =  /  F.dv,  (00) 

Jo  Jm 

where  v0  is  the  low-frequency  limit  for  photo- 
electron  pnxhiction.  Equation  (00)  has  lieen 
evaluated  in  three  parts,  corresponding  to 
regions  effective  in  ionization  of  X  +  X2, 
O,  and  ()•„>,  respectively,  since  the  great 
difference  in  altitude  of  ion  formation  leads 
one  to  expect  considerably  different  values 
of  the  recombination  coefficient  for  each 
region.  Substituting  numerical  values  in 
order  gives  \qP{h)  dh  =  [1.8  X  10“  +  1.4  X 
It)11  +  2  X  10"J  cm-2  sec-1. 

The  formation  of  an  ionosphere  is  con¬ 
siderably  more  complicated  than  presented 
here.  It  should  be  emphasized  that  the 
quantities  derived  in  this  section  are  correct 
only  to  order  of  magnitude.  It  is,  however,  of 
interest  to  note  that  the  derived  quantities 
are  of  the  same  order  as  those  due  to  the  solar 
proton  flux. 


The  Effective  Recombination 
Coefficient 

In  general,  electrons  and  ions  can  recom¬ 
bine  by  the  following  processes: 

Radiative  recombination  (atomic  ions) 

X+  +  e  —  X'  +  hv  (m) 

Dielectronie  recombination  (atomic  ions) 
X+  +  c  — X'— *X'  +  hv  (n) 

Dissociative  reeombinat  ion  (molecular 
ions) 

XV+  +  e  — >  X'  +  Y'  (p) 

Three-body  recombination 

X+  +  e  +  Y  — ►  X  +  Y  (r) 

or 

XY+  +  e  +  /  —  XY  +  Z,  (s) 

where  the  prime  denotes  that  the  atom  or 
molecule  is  in  an  excited  state.  The  double 
prime  indicates  an  excited  state  lying  in  the 
continuum.  A  discussion  of  these  processes 
is  given  by  Xawrocki  and  Papa  (1901), 
whose  tabulated  rate  coefficients  are  adopted 
in  the  following  discussion. 

Radiative  recombination  coefficients  ar 
measured  for  atomic  oxygen  vary  from  3.4  X 
10-12  cm3  sec-1  at  T  =  230°K  to  0.8  X  10-12 
cm3  sec-1  at  T  =  2000°K.  Since  the  major 
contribution  to  the  coefficient  is  due  to 
captures  into  excited  states,  it  is  likely  that 
similar  values  will  apply  to  nitrogen  radia¬ 
tive  recombination. 

Rate  coefficients  for  dielectronie  recom¬ 
bination  are  probably  on  the  order  of  10-10 
to  10-12  cm3  sec-1;  however,  only  meager 
data  is  available  on  this  process. 

Dissociative  recombination  coefficients  a,i 
for  X2+  have  measured  values  which  range 
from  1.4  X  10-7  cm3  sec-1  to  8.5  X  10-7 
cm5  sec-1  at  room  temperature,  and  1.4  X 
10-6  cm3  sec-1  at  T  =  2500°K.  It  seems  that 
for  the  present  discussion  a  value  of  10-7 
cm3  sec-1  is  a  conservative  choice.  Values  for 
02+,  XO+,  C()2+,  and  II20+  are  all  about 
10-7  cm3  sec-1  with  an  uncertainty  of  ±2 
in  the  exponent. 
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Three-body  rat(‘  coefficients,  0,  have  typi¬ 
cal  values  of  5  X  10  ■'  cm*  sec-1  for  X2+  and 
()2+  and  10-si±I  cm®  sec-1  for  other  reactions 
involving  molecules.  Similar  values  an?  at¬ 
tained  for  three-lnxlv  recombination  of 
atomic  ions. 

From  the  data  of  Tables  IN'  and  Y  it  is 
seen  that  the  maximum  depth  of  proton 
penetration  is  0.01  cm-atm  for  the  models 
considered.  This  depth  thus  forms  the  lower 
boundary  to  the  proton  ionization  layer. 
At  this  level  Eq.  (20)  becomes 

Hn(hmi„)  =  2.8  X  10l2m/7',  (07) 


which  for  the  range  of  m/T  considered  in 
the  models  yields  a  neutral  particle  num¬ 
ber  density  at  the  lower  limit  of  proton 
ionization  in  the  range  2  X  10"  cm-3  to 
8  X  10,#  cm-3. 

With  the  aid  of  Eq.  (34)  we  may  compare 
the  three-body  recombination  rate  to  that 
of  radiative  recombination.  Three-body  re¬ 
combinations  will  dominate  if,  for  n,-  =  n„ 
the  total  particle  density  N  »  ar/0.  From 
the  above  discussion  ar/0  ~  1015  cm-3;  thus 
for  three-body  reactions  to  be  important  N 
must  be  »  1015  cm-3.  It  should  be  recalled 
that  N  at  the  occupation  level  is  about  1014 
cm3.  One  may  conclude  then  that  three-body 
electron-ion  recombination  processes  are  not 
effective  above  the  occultation  level. 

The  rate  of  radiative  recombination  de¬ 
pends  upon  the  number  density  of  atoms 
n( A),  while  that  for  dissociative  recombina¬ 
tion  on  the  number  density  of  molecules 
n (M).  If  we  neglect  ion-atom  interchange, 
we  may  compare  these  rates  by 


Radiative  recombination  rate 
Dissociative  recombination  rate 

-  a'nWn' 

a„n(  M)«, 


(68) 


These  rates  are  equal  for  n(A)/«(M)  = 
a,//ar  ~  105.  It  is  apparent  that  a  very  high 
degree  of  dissociation  is  required  for  radia¬ 
tive  recombination  to  be  important.  A 
similar  conclusion  may  be  reached  concern¬ 
ing  dielcctronic  recombination,  although  if 
the  rate  coefficient  can  be  as  large  as  10-,° 
cm3  sec-1  the  rates  of  dielectronic  and  dis¬ 
sociative  recombination  would  be  equal  for 


n(A)/n(M)  ~  100,  a  much  smaller  degree  of 
dissociation. 

To  achieve  this  high  degree  of  dissociation 
with  the  assumed  large  abundance  of  X* 
requires  that  the  Herzberg-Herzberg  mecha¬ 
nism  lx?  very  efficient.  This  requires  a  low 
abundance  of  H,  H2,  02,  Oj,  and  XO  since 
these  molecules  liave  moderate  to  strong 
absorptions  in  the  region  1 200-1 250A. .  Colli- 
sional  frequencies  of  molecules  and  atoms 
in  the  region  of  proton  ionization  are  on  the 
order  of  1  per  sec;  thus  collisional  de-exeita- 
tion  of  X»  in  the  “Herzberg  State”  prior  to 
dissociation  is  unlikely  (lifetime  of  state  ~ 
10~3  sec). 

The  dissociation  of  X2  in  the  terrestrial 
atmosphere  through  the  Herzberg-Herzberg 
process  is  limited  by  diffusion  (see,  c.g., 
Bates,  1954),  and  it  is  easy  to  show  that  a 
similar  conclusion  applies  to  Venus.  A 
simple  solution  of  the  one-dimensional  dif¬ 
fusion  equation  yields  for  the  diffusion 
current  of  the  yth  molecular  species, 

J  i  —  —  Dj(AHj/Ah),  (69) 

where  AnffAh  is  the  concentration  gradient, 
and 


^  &kT/m,y» 

12t  nr-  ' 


(70) 


A  characteristic  time  for  the  establishment 
of  diffusive  equilibrium  in  an  isothermal 
region  above  a  critical  level,  hc,  characterized 
by  an  overall  number  density  nc,  and  a 
number  density  nCJ  for  the  jlh  constituent  is 

t,  =  nciH,/Jh  (71) 


where  Hj  is  the  scale  height  of  the  jth  con¬ 
stituent  in  the  diffusion  regime.  Substitution 
of  (69)  and  (70)  into  (71)  yields  for  the 
characteristic  velocity  of  the  jth  constituent, 
diffusing  between  some  altitude  h  =  Ah  +  hc 
where  n  « nc  and  the  critical  level,  the 
expression 


_ V3/jlYb 

e]  tj  —  I2imcr*  \HjJ  ' 


(72) 


Diffusive  equilibrium  exists  when  vCj  is  much 
larger  than  the  prevailing  turbulent  veloci¬ 
ties.  For  the  Earth,  diffusive  equilibrium  is 
established  when  vej-  exceeds  a  few  cm  sec-1. 


120 


It.  CJ.  WALKEK  AND  C.  SAGAN 


The  factors  controlling  the  prevailing  turbu¬ 
lent  velocities  in  the  Earth's  upper  atmos¬ 
phere  are  very  complicated  and  almost 
entirely  unknown.  Venus  is  a  planet  which 
is  very  slowly  rotating.  Some  very  simple 
models  of  the  atmospheric  circulat  ion  for  the 
likely  high  surface  pressures  give  prevailing 
wind  velocities  which  arc  not  excessive  by 
terrestrial  standards  (Mintz,  19<>2).  In  the 
absence  of  any  further  information  we  will 
adopt  the  same  value  for  *v>  for  the  estalx- 
lishment  of  diffusive  equilibrium  that  pre¬ 
vails  on  Earth.  The  qualitative  conclusions 
which  we  will  later  draw  in  this  paper  do  not 
depend  critically  on  this  assumption.  We 
then  find, 

itr  3  X  10“  cm-* 

for  II ,  ~  30  km.  This  value  is  not  sensitively 
dependent  on  the  choice  of  Hj.  Thus,  dif¬ 
fusive  equilibrium  is  established  at  about 
the  level  of  proton  and  ultraviolet  ionization, 
but  significantly  above  the  occultation  level. 

The  X2  dissociation  rate  through  the 
Herzlierg-IIerzlierg  mechanism  on  Venus  is 

q( j)  ~  4  X  10-,i«(X2). 

Taking  the  atmosphere  to  be  primarily  X2 
at  the  ionization  level,  we  find  then* 

</(j)  ^  0.4  cm-*  sec-'. 

The  mean  dissociation  time  is  then  10"/0.4 
~  2  X  10"  sec,  or  several  thousand  years. 
Iiut  the  corresponding  times  for  diffusive 
mixing  of  the  atomic  nitrogen  photodissocia¬ 
tion  products  of  the  Herzberg-1  lerzberg 
mechanism  is,  by  (71),  of  the  order  of  days. 
When  the  nit  rogen  atoms  diffuse  to  altitudes 
when*  the  total  number  density  is  10"  or 
1013  cm-3,  they  rapidly  recombine  (s  >e,  e.g., 
Bates,  1054)  by  the  three-body  process  (1). 
Accordingly,  the  Herzlx>rg-Hcrzlx*rg  mecha¬ 
nism  is  generally  inefficient  on  Venus; 
molecular  nitrogen  will  lx*  abundant  at  the 
ionization  level;  and  dissociative  recombina¬ 
tion  cannot  be  avoided. 

We  conclude  that  the  effective  recombina¬ 
tion  coefficient  near  the  region  of  maximum 
ionization  in  the  upper  Cytherean  atmos¬ 
phere  is 

ae  ~  10-7  cm*  sec"  ’.  (73) 


Calculation  of  the  ionization  rate  by 
protons  involved  the  assumption  that  the 
electron  temperature  was  equal  to  the  local 
kinetic  temperature.  The  validity  of  this 
assumption  will  now  lx;  examined.  The  time 
/  for  electrons  eject cd  with  excess  energy  I ’ 
(electron  volts)  to  thermalize  at  temperature 
T  is  given  by  Ginzburg  (1901)  as 


/  = 


In  (r//.- 7’) 


5 


sec, 


(74) 


where  f:  is  Boltzmann’s  constant  and  for 
clast  ic  collisions  in  atomic  gases  5  =  1  .OSS  X 
10-*  .1,  .1  lx*ing  the  atomic  weight  of  the 
particles.  For  inelastic  collisions  in  air  and 
similar  molecular  mixtures  5  ~  10~*  [see 
Ginzburg  (1901),  Table  5.1.  p.  09  for  exact 
values]. 

From  Eq.  (29)  the  average  lifetime  of  an 
ion  is  At  =  (Ar,ac)~]  sec,  where  AT,-  is  the 
number  density  of  ions.  If  At  >  t  the  elec¬ 
trons  will  have  time  to  thermalize  and 
Te  =  7’;  thus  for  electrons  to  thermalize  we 


requin* 


ar  < 


v,b 

cTTh TtJT/W)' 


(75) 


It  can  be  easily  shown  that  for  ft*  <  2  X  106 
E<|.  (22)  determines  vr.  Substitution  of 
(22),  5  =  10-*,  and  V  =  A’,  —  Ionization 
potential  ~  15  eV  yields  the  requirement 
that  a,  <2X  10-7  cm*  sec-1.  Since  tin* 
deduced  value  of  ae  is  10-7  cm*  sec-1  tin* 
assumption  Te  =  T  is  approximately  valid. 

Results  and  Conclusions 

From  t he  results  displayed  in  Tables  V 
and  VI,  it  is  seen  that  the  maximum  rates 
of  ionization  consistent  with  the  observed 
proton  spectrum  an*  X  10*  cm-*  sec-1, 
while  the  upper  values  of  Z  are  ~5  X  106 
cm-3  sec-1.  From  Kqs.  (51)  and  (53),  we  find 
jne-dh  ^7  X  1017  cm-5  for  a,  ~  10-7  cm* 
sec-1.  In  all  our  proton  ionization  models, 
the  ionized  layers  are  ~5t)  km  thick  'of. 
Figs.  1-4).  The  corresponding  mean  elec,  on 
density  is  ~  4  X  10s  cm-*.  Phot oionizat  ion 
will  contribute  an  electron  density  of  the 
same  order  of  magnitude.  Had  radiative 
recombination  and  a, .  ~  10-13  cm’  sec-1 
been  realizable,  values  of  fnc-  dh  approaching 
1023  cm-5  and  values  of  n,  ~  10*  cm-*  would 
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have  })een  derived.  Even  these  electron 
densities  are  too  low  for  the  ionospheric 
model.  Hut  liecause  of  the  great  abundance 
of  X2  in  the  region  where  ionization  occurs, 
values  of  ne  even  as  large  as  107  cm-3  are 
excluded. 

For  Z  =  5  X  10®  cm-2  sec-1  and  a,  =  10-7 
cm3  sec-1,  the  optical  depth  of  the  ionized 
layer,  due  to  free-free  emission,  would  be, 
from  Eqs.  (9)  and  (5(5),  r(X)  ~  4  X  10-9  X2. 
Significant  opacity  occurs  only  for  X  >  100 
meter,  far  beyond  the  microwave  region  of 
relevance  for  Venus.  Had  we  chosen  ar  ~ 
10-12  cm3  sec-1,  we  would  have  found  r(X)  ~ 
4  X  10-4  X2,  implying  t(10  cm)  =  0.04,  and 
r((58  cm)  =  1.8.  Even  these  optical  depths 
are  inadequate  to  explain  the  microwave 
observations  by  free-free  emission,  although 
they  would  have  caused  significant  attenua¬ 
tion  of  long-wavelength  radar  returns  from 
Venus. 

The  models  developed  for  the  formation 
of  a  Cytherean  ionosphere  by  solar  protons 
assumed  a  vanishing  magnetic  field  strength. 
It  is  well  known  that  interactions  between 
a  neutral  plasma  and  an  external  magnetic 
field  are  governed  by  the  relative  energy 
densities  of  the  field  and  the  plasma.  If  the 
energy  density  of  the  plasma  is  greater  than 
that  of  the  field,  charge  separation  does  not 
occur  and  the  trajectories  of  the  individual 
particles  are  governed  by  the  plasma  rather 
than  by  the  magnetic  field.  The  energy 
density  Up  of  the  proton  stream  is 

Up  =  \m  ^  uojVj2  =  V  iiajKj.  (70) 

i  J 

The  energy  density  Um  of  the  magnetic  field 
is 

Un  =  (77) 

If  we  require  Cp  >  Um  Hie  proton  stream 
will  not  be  appreciably  deflected  by  (he 
magnetic  field;  this  condition  can  be  written 

B  <  (8tt  ^  no MY12-  (78) 

> 

Substitution  of  numerical  values  from  Table 
IV  yields  B  <  1.04  X  10-3  gauss. 

The  Mariner  2  magnetometer  experiment 
(Smith  d  ah,  1903)  detected  no  sign  of  a 


Cytherean  magnetic  field  at  a  distance  of 
4.1  X  104  km  from  the  planet.  This  result, 
and  similar  negative  results  in  the  search  for 
Cytherean  radiation  belts,  lead  to  the  con¬ 
clusion  that  the  magnetic  moment  of  Venus 
is  less  tlian  that  of  Earth;  however,  a  surface 
magnetic  field  strength  of  0.1  gaass  is  not 
inconsistent  with  the  observations  [see,  e.g., 
Sonett  (1903)].  The  slow  rotation  rate  of 
Venus,  obtained  from  radar  observations,  is 
also  consistent  with  a  low  magnetic  field, 
but  theories  of  the  origin  of  the  geomagnetic 
field  are  not  nearly  in  a  sufficiently  satisfac¬ 
tory  state  to  predict  whether  the  Cytherean 
surface  field  strength  is  greater  than  or  less 
than  10-3  gauss. 

The  effect  of  a  local  magnetic  field  on  the 
microwave  absorption  coefficient  and  thus 
on  the  optical  depth  of  the  ionospheric  layer 
is  well  known.  Magnetic  field  effects  become 
important  when  the  gyrofrequency  ap¬ 
proaches  the  angular  frequency  w  of  the 
wave;  then  [see,  e.g.,  Ratcliffe  (1059)] 
~  10-4  B\,  when  X  is  in  cm  and  B  in 
gauss.  For  X  =  (58  cm ,un/u  =  1  for  B  —  147 
gauss,  a  value  not  admitted  by  the  Mariner 
2  observations. 

Thus,  no  permissible  combination  of 
Cytherean  magnetic  field  strength  and 
atmospheric  composition  yields  results  com¬ 
patible  with  the  ionospheric  model.  If  we 
refuse  to  invoke  some  other  unknown  ioniza¬ 
tion  source,  it  follows  that  we  must  reject 
the  ionospheric  model.  A  similar  conclusion 
was  reached  by  Sagan,  Siegel,  and  Jones 
(19(51)  and  by  Kellogg  and  Sagan  (19(51). 
Our  maximum  derived  electron  densities 
are  much  smaller  than  those  required  by 
Jones  (19(51),  by  Priester,  Roomer,  and 
Schmidt-Kaler  (19(52),  and  by  Kuzmin 
(19(54).  The  discussion  of  Danilov  and 
Yatscnko  (19(53)  requires  that  the  iono¬ 
spheric  optical  depth  at  10  cm  be  of  the 
order  of  2,  a  value  much  larger  than  the 
known  ionization  sources  permit. 

If  we  were  to  postulate  additional  proton 
fluxes  with  energies  outside  the  range  meas¬ 
ured  by  Mariner  2,  we  would  be  laced  with 
the  following  difficulties: 

1.  If  an  additional  flux  of  high-energy 
protons  is  postulated,  these  protons  will 
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penetrate  deeply  into  the  atmosphere  and 
produce  most  of  their  ionization  in  regions 
where  the  recombination  coefficient  is  large, 
thus  greatly  reducing  their  effectiveness. 

2.  If  a  flux  of  low-energy  particles  is  postu¬ 
lated,  the  required  flux  must  be  increased 
many  orders  of  magnitude  to  make  up  for  the 
loss  in  ionization  efficiency  of  the  low-energy 
protons. 

The  question  of  the  proton  flux  at  smaller 
energies  than  those  measured  by  Mariner  2 
is  relevant  to  the  observed  correlation  of  the 
68-cm  radar  determination  of  the  distance 
to  Venus  with  the  10.7-cm  and  20-cm  solar 
flux.  No  such  correlation  has  been  reported 
for  the  12.5-cm  radar  ranging.  Muhleman 
(1963)  has  postulated  that  solar  activity 
somehow  reduces  the  Cytherean  electron 
density,  temporarily  increasing  the  group 
velocity  of  the  radar  pulse  in  the  Cytherean 
ionosphere,  and  decreasing  the  measured 
value  of  the  astronomical  unit.  Muhleman 
can  match  the  observations  if  fne2  dz  ~  5  X 
1022  cm-5,  provided  nt  ~  107  cm-3  and  Az  ~ 
104  km.  This  assumed  ionospheric  scale 
height  is  much  greater  than  the  measured 
solar  proton  flux  admits  (ef.  Figs.  1-4,  where 
A2~50  km).  Muhleman ’s  explanation  of 
the  correlation  between  the  radar  range  and 
solar  activity  is,  therefore,  admissible  only 
if  there  exists  an  alternative  ionization 
source — perhaps  a  very  large  flux  of  low- 
energy  protons  which  are  thermalized  at 
very  high  altitudes  where  radiative  recom¬ 
bination  (due  to  0+,  for  example)  dominates. 
But  even  then,  the  admissible  values  of 
fne2  dz  are  inconsistent  with  the  ionospheric 
model  of  the  microwave  emission,  as  Muhle¬ 
man  (1903)  himself  emphasized. 

The  foregoing  difficulties  with  the  iono¬ 
spheric  model  lead  us  to  conclude  that  the 
surface  of  Venus  is  the  source  of  the  observed 
microwave  emission.  Greenhouse  models  of 
Venus  now  appear  capable  of  explaining 
the  high  surface  temperatures  and  the  ob¬ 
served  microwave  spectrum  and  phase  effect 
(Sagan  and  Pollack,  196.5;  Pollack  and 
Sagan,  1966). 

Acknowledgment 

We  are  grateful  to  Prof.  David  R.  Bates  for  a 
critical  reading  of  the  manuscript. 


Refeuexces 

Baratii,  F.  T.,  Barrett,  A.  H.,  Copeland,  J., 
Jones,  D.  E.,  and  Lilley,  A.  E.  (1963).  Mariner 
II:  Preliminary  reports  on  measurements  of 
Venus — Microwave  radiometers.  Science  139, 
90S-909. 

Barrett,  A.  H.,  and  Staelin,  D.  H.  (1964).  Radio 
observations  of  Venus  and  the  interpretations. 
Space  Sri.  Revs.  3,  109-135. 

Bates,  D.  R.  (1954).  The  physics  of  the  upper 
atmosphere.  In  “The  Earth  as  a  Planet”  (G.  P. 
Kuiper,  ed.),  pp.  576-636.  Univ.  of  Chicago 
Press,  Chicago,  Illinois. 

Bottema,  M.,  Plummer,  W.,  and  Strong,  J. 
(1964).  Water  vapor  in  the  atmosphere  of 
Venus.  Astrophys.  J.  139,  1021-1022. 

BranuT,  J.  C.  (1962).  A  model  of  the  inter¬ 
planetary  medium.  Icarus  1,  1-6. 

Chamberlain,  J.  W.  (1965).  Astrophys.  J.  141, 
1184-1205. 

Cook,  C.  J.,  Jones,  E.,  Jr.,  and  Jorgensen,  T., 
Jr.  (1953).  Range-energy  relations  of  10-  to 
250-kcV  protons  and  helium  ions  in  various 
gases.  Rhys.  Rev.  91,  1417-1422. 

Danilov,  A.  D.  (1964).  Radioastronornical  re¬ 
search  and  current  concepts  of  the  atmosphere 
of  Venus.  Kosmichoskiyr  Isslcdovaniya  2,  188- 
215. 

Danilov,  A.  D.,  and  Yatsenko,  S.  P.  (1963).  An 
ionospheric  interpretation  of  results  of  radio 
observations  of  Venus.  Parts  I  and  II.  Geo¬ 
magnetism  and  Acronomy  3,  585-597. 

de  Vaucoui.eurs,  G.,  and  Menzel,  D.  H.  (1960). 
Results  of  the  occultation  of  rcgulus  by  Venus. 
Nature  188,  28-33. 

Drake,  F.  D.  (1962).  10-cm  observations  of  Venus 
in  1961.  Publ.  Nall.  Radio  Astron.  Obs.  1, 
165-178. 

Ginzburg,  V.  L.  (1961).  “Propagation  of  Electro¬ 
magnetic  Waves  in  Plasma”  (W.  L.  Sadowski 
and  D.  M.  Gallik,  eds.).  Gordon  and  Breach, 
New  York. 

Jones,  D.  E.  (1961).  The  microwave  temperature 
of  Venus.  Planetary  Space  Sci.  5,  166-167. 

Kellogg,  W.  W.,  and  Sagan,  C.  (1961).  The 
atmospheres  of  Mars  and  Venus.  Natl.  Acad. 
Sci-Natl.  Res.  Council  Publ.  944  (Washington, 
D.  C.). 

Kotelnikov,  V.  A.  (1961).  Radar  contact  with 
Venus.  J.  Brit.  Inst.  Radio  Engrs.  621,  293-295. 

KouRGANOrF,  V.  (1952).  “Basic  Methods  in 
Transfer  Problems.”  Clarendon  Press,  Oxford, 
England. 

Kuzmin,  A.  D.  (1964).  On  the  ionospheric  model 
of  Venus.  In  “Life  Sciences  and  Space  Re¬ 
search,  II”  (M.  Florkin  and  A.  Dollfus,  eds.), 


IONOSPHERIC  MODEL  OF  THE  VENUS  MICROWAVE  EMISSION 


123 


pp.  211-221.  North-Holland  Publ.  Co.,  Amster¬ 
dam. 

Mayer,  C.  H.  (1903).  General  report  on  planetary 
radio  astronomy.  Proc.  Intern.  Aslrophys.  Col- 
loq.,  Illh,  Liege,  1963,  pp.  99-111. 

Mintz,  Y.  (1962).  The  energy  budget  and  atmos¬ 
pheric  circulation  on  a  synchronously  rotating 
planet.  Icarus  2,  172-173. 

Minzner,  R.  A.,  and  Ripley,  W.  S.  (1957). 
"Handbook  of  Geophysics,”  1st  ed.  Air  Force 
Cambridge  Res.  Lab.,  Bedford,  Massachusetts. 

Mitka,  S.  K.  (1952).  “The  Upper  Atmosphere,” 
2nd  ed.  The  Asiatic  Society,  Calcutta,  India. 

Muiii.eai.an,  D.  O.  (1963).  The  electrical  charac¬ 
teristics  of  the  atmosphere  and  surface  of  Venus 
from  radar  observations.  Icarus  1,  401-411. 

Xawrocki,  P.  J.,  and  Papa,  R.  (1961).  Atmos¬ 
pheric  processes.  Air  Force  Cambridge  lies. 
Lab.  Rcpl.  No.  AFCRL-595. 

Xeugeuauer,  M.,  and  Snyder,  C.  W.  (1962).  Pre¬ 
liminary  results  from  the  Mariner  2  solar 
plasma  experiment.  Jet  Propulsion  Lab.,  Calif. 
Inst.  Tcchnol.,  Peek.  Memo.  No.  33-111. 

i)i‘iK,  E.  J.  (1961).  The  aeolosphere  and  atmos¬ 
phere  of  Venus.  J.  Gcophys.  lies.  66,  2807- 
2819. 

Oster,  L.  (1961).  Emission,  absorpti  m,  and  con¬ 
ductivity  of  a  fully  ionized  gas  at  radio  fre¬ 
quencies.  Rev.  Mod.  Phys.  33,  525-5-13. 

Pettf.ngiu.,  G.  H.,  Briscoe,  H.  W.,  Evans,  J.  V., 
Gkhrels,  E.,  Hyde,  G.  M.,  Kraft.  L.  G.,  Price, 
R.,  and  Smith,  W.  B.  (1962).  A  radar  inves¬ 
tigation  of  Venus.  Aslron.  J.  67,  181-190. 

Pollack,  J.  B.,  and  Sagan,  C.  (1965).  The  inicro- 


Avave  phase  effect  of  Venus.  Icarus  4,  62-103. 

Priester,  W,  Roemer,  M.,  and  Schmidt-Kaler, 
Tii.  (1962).  Apparent  relation  between  solar 
activity  and  the  440-Mc/s  radar  distance  of 
Venus.  Nature  196,  464-465. 

Ratcliffe,  J.  A.  (1959).  “Magnetoionic  Theory.” 
Cambridge  Univ.  Press,  Cambridge,  England. 

Sagan,  C.  (1960).  The  radiation  balance  of  Venus. 
Jet  Propulsion  Lab.,  Calif.  Inst.  Technol.,  Tech. 
Kept.  32-34. 

Sagan,  C.  (1961).  The  planet  Venus.  Science  133, 
849-858. 

Sagan,  C.  (1962).  Structure  of  the  lower  atmos¬ 
phere  of  Venus.  Icarus  1,  151-169. 

Sagan,  C.,  and  Kellogg,  W.  W.  (1963).  The  ter¬ 
restrial  planets.  Ann.  Rev.  Aslron.  Aslrophys. 
1,  235-266. 

Sagan,  C.,  and  Pollack,  J.  B.  (1965).  On  the 
nature  of  the  clouds  and  the  origin  of  the 
surface  temperature  of  Venus.  To  be  published. 

Sagan,  C.,  Siegel,  K.  M.,  and  Jones,  D.  E.  (1961). 
On  the  origin  of  the  Venus  microwave  emis¬ 
sion  (abstract).  Aslron.  J.  66,  52. 

Smith,  E.  J.,  Davis,  L.,  Coleman,  P.  J.,  and 
Sonett,  C.  P.  (1963).  Mariner  2:  Preliminary 
measurements  of  Venus  —  Magnetic  field. 
Science  139,  909-910. 

Sonett,  C.  P.  (1963).  A  summary  review  of  the 
scientific  findings  of  the  Mariner  Venus  mis¬ 
sion.  Space  Sci.  Revs.  2,  751-777. 

Victor,  W.  K.,  and  Stevens,  R.  (1961).  Explora¬ 
tion  of  Venus  by  radar.  Science  134,  46-48. 

Whipple,  F.  L.  (1962).  Quoted  by  Muhleinnn 
(1963). 


